HU  HD. 

DDC  FILEJQea  AD  AO  58086 


Finite  Element  Models  with  Velocity  Discontinuities 


Hendrik  M.  van  Rij,  Research  Assistant 

and 

Philip  G.  Hodge  Jr.,  Professor  of  Mechanics 


Department  of  Aerospace  Engineering  and  Mechanics 
University  of  Minnesota 
Minneapolis,  Minnesota  55455 


/ 


n n 

> iJ 


nxsaxsmaiTJFS  ksst  quality  p -AonoMAa 

Vto  COPY  TURWISH©  TO  EDC  OWTAIHED  A / 

r'mmcMz  maim  o?  faobs  rarcHDO*o* 

s®W>EWCS 


■ 

! AUG  25 


May  1978 


; . : ici  ’asei 


Technical  Report 


/ 


I 


Prepared  for 


OFFICE  OF  NAVAL  RESEARCH 
Arlington,  VA  22217 

OFFICE  OF  NAVAL  RESEARCH 
Chicago  Branch  Office 
536  South  Clark  3t. 
Chicago,  IL  60605 


78  07  12  065 


H 


DISCLAIMER  NOTICE 


THIS  DOCUMENT  IS  BEST  QUALITY 
PRACTICABLE.  THE  COPY  FURNISHED 
TO  DDC  CONTAINED  A SIGNIFICANT 
NUMBER  OF  PAGES  WHICH  DO  NOT 
REPRODUCE  LEGIBLY. 


JECUHITY  CLAUDICATION  or  THIS  »»Ct  Of  rjy <•*•</) 


READ  TR3TRUCTIONS 
BEFORE  CO'JPLETTNC.  FORM 


2.  GOVT  ACCESSION  KOJ  1-  RECIPIENT**  CATALOG  NUMBER 


REPORT  CCCIWENTAT1CN  PAGE 


■.  TY*e  of  report  * rtwoo  covzrco 


STRUCTURAL  ^ELASTICITY  Xg 

Finite  Element  Models  with 
Velocity  Discontinuities* 


m 


P r of  "f  I1  11  frCS~~ 


IMp  if— 

n5t-75-C-J0177j 


UUM«ERf»> 


t.  PERFORMING  ORGANIZATION  NAME  AND  ADORIJJ 

University  of  Minnesota 
Minneapolis,  Minnesota-  55455 


II.  CONTROLLING  OFFICE  NAME  AMD  ADDRE4S  ... 

OFFICE  OF  NAVAL  RESEARCH  . .. 

Chicago  Branch  Office 

___L36^uth£!££k_JLt*— z— Clli£IlQ£l*— 1L— 6H£H5_ 

14.  MONITORING  AGENCY  NAME  A ADORES*/!  dlllonnl  tram  Conlrolllnt  Ollleo ) 


10.  PROGRAM  ELEMENT.  PROJECT,  TAJK 
AREA  A BORA  UNIT  NUMBERS 


£ NR  064-429 


K3H 


IS.  NUMB 


[fCLloS,  (at  th!a  rtpari) 


Unclassified 


It.  DISTRIBUTION  STATEMENT  (ol  thlo  Report) 


17.  DISTRIBUTION  STATEMENT  «A»«r*«l  onlond  In  Blook  30,  II  dlllotonl  tram  Report)  J- 


IS.  supplementary  notes 


IB.  KEY  PONDS  (Conllnvo  on  nwn  ode  It  roeoooory  tnd  Idonlll y Of  block  nmmbof) 


Plasticity,  perfectly-plastic , finite-elements,  plane  strain, 
Prandtl  punch,  slip,  velocity  discontinuities,  notched  bar.  . 


so.  abstract  (Continue  an  .<  <«  a n»o*f~T  and  Identify  by  noo*  nvmtoff  Conventional  finite-element 

models  are  based  on  displacement  or  velocity  fields  which  are  at  least  continuous 
However,  it  is  known  that  perfectly  plastic  materials  may  involve  discontinuities 
of  the  tangential  velocity  component  along  certain  lines.  In  this  investigation, 
two  two-dimensional  finite  element  models  are. proposed  which  will  allow  for 
such  discontinuities.  • . • . . .*  * . • 

A regular  pattern  of  i^iangular  elements  is  assumed.  In  the  first  model,  the 
elements  are  assumed/ to\be  rigid,  and  in  the  second  model,  we  have  a linear  , 


displacement  field  in  each  element.  Across  the  line  separating  any  two 
adjoining  elements  the  normal  displacement  component  is  continuous,  but 
a discontinuity  may  exist  in  the  tangential  component.  The  defining 
equations  - compatibility,  equilibrium,  and  constitutive  - are  developed 
with  the  aid  of  the  Principle  of  Virtual  Work. 

Prandtl's  punch  problem  and  tension  in  notched  bars  are  solved  under 
plane  strain  conditions.  Comparison  is  made  with  existing  analytical 
and  other  numerical  solutions,  in  order  to  evaluate  the  merits  of  allowing 
for  discontinuities. 


\ 


A general  approach  to  any  finite  element  model 


inbers  in  square  brackets  refer  to  items  listed  in  Hodge  (£]  showed  that  when  this  method  is  applied 

is  bibliography. 

to  the  constant-strain  triangle  model  for  two-dimensional 


first  obtained  by  Turner 


\ 


model  which  will  allow  for  such  discontinuities 


rig.  2 and  d the  displacement  of  the  triangle  AOO  along 


involved  with  different  small  node  mechanisms,  so  placement  of  the  triangle  above  PQ.  Then 


Ill]  for  generalised  stress  and  strain  IF  |<J  | < k THEN  & “ G t (2-lla) 


Eq.  (2-12)  as  will  be  illustrated  in  Chapter  V. 


alternative  viewpoint,  whicn  is  more  useful 


examples  to  be  considered  later,  all  domains 


Note  that  if  the  edge  were  longer  thi  next  node  would  where  Ti  is  the  averaged  dimensionless  vertical  force 


equilibrium  equation  (2-16a)  and  (2-18)  then  serve  merely  as  definitions 


checked  aqainst  the  requirement  that  w-_  £ 0 and 


The  stress  rates  are  all  uniquely  determined 


1 


of  tha  output,  rather  than  by  the 


possible  physical  significance  of  this  alternating 


yield.  Finally,  it  yields  a third  time  to  terminate  causes  a zero  in  the  denominator  of  Eg.  <2-28b) , but 


The  agreement  of  both  numerical  solutions  with 


(S-E) 


solution  (14,15)  for  the  rough  punch  on  an  infinite 


stopped  and  the  combined  mechanism  shown  in  Fig. 


r 


<0 


and  in  all  cases  different  from  the  Prandtl  punch  Classical  Model 

problem,  even  though  the  present  problem  is  frequently  Traditionally  finite  element  models  are  based  or. 

viewed  as  the  tensile  analogue  of  the  Prandtl  punch.  on  nmnfn  u<>h<n 


consistent  field  of  a rigid  perfectly-plastic  material 


obtained  in  similar  fashion.  The  relations 


Tig.  38 


Eqa.  (5-10)  can  be  solved  for  the  stress  rates 


Continuation  of  this  procedure  leads  to  Eqs.  (5-18)  were  derived  from  Eq.  (5-14)  by  using 

_ the  slip  relations.  However,  they  could  also  be 

»1  - U - l vB  ♦ *0(V8H,1V(*/I) 

written  down  by  symmetry  considerations.  It  follows 


the  kinematics  of  the  combined  model  are  fully 


If  <JB„  denotes  the  boundary  slip,  then 


the  boundary  forces  T,  and  or  in  terms  of  the 


velocity  discontinuities  are  not  allowed 


behavior  of  the  triangular  elements  in  that  more  than 


at  the  yield  point  load.  In  view  of  these 


the  load  deflection  curve  for  the  classical  model  up 


punch,  and  Figs.  <2  and  46  are  very  much  alike 


elastically  undetermined  unit  mechanism  at  and  the  double  edge  notched  problem  is  that  the 


the  same  fashion  as  we  shaded  Fig.  44  for  the  rough  yielding  mechanism  at  which  the 


As  with  the  double-notched  bar,  there  is  only 


Chapter  VIII  Conclusions  I that  of  an  analytic  solution  for  a semi -infinite 

domain  with  v « 0.50  [14,15]  shown  in  Fig.  15  except 


lassical  models  are  the  Eame.  The  displacement  pattern 


the  sexi-infinite  domain,  Prandtl  ( 9 ) obtained  in  Fig.  17. 


product  with  the  vector  of  basic  unknowns  is  set  equal  rate  for  each  yielding  element  is  positive.  If  it  is 


Based  upon  very  limited  information,  the  slip  model  bY  deleting  the  slip  part.  This  is  due  to  the  fact 

appears  to  be  the  most  efficient  in  terms  of  CPU  time.  that  the  kinematic  variables  associated  with  the  slip 


extracted  in  this  range  from  the  combined  model  results  , made  on  unloading.  The  answers  are  acceptable  up  tc 


\ 


edge  expansion,  the  singularities  in  the  stiffness 


Finite  Element  Method",  van  Nostrand 


8th  Intern.  Congr.  Appl.  Kech.",  Istanbul 


100 


Using  the  geometrical  relations  Substituting  45*  in  all  angles  of  Eq.  (A-l)  we 


Z01  - " T0T 


I 


I 


103 


which  the  subroutines  SLIP  and  REG  are  written.  As  part  the  adjustment  in  the  stiffness 


w 


) 23  — / */ 


instability  without  reaching  a limit  load,  it  seemed 


to  taking  a Poisson  ratio  of  .49999  instead  of 


wo  dropped  the  superscripts  as  element  identifi- 


incorporated  in  the  stiffness  matrix.  As  convergence  PPFLOrfciNPcTtTAPC^iNPuT.TAPEft.ouiPuwccatUPFi.uPtii 


.232  is  reliable  for  the  Prandtl  rough  punch  with 


ITERATION  NU*BCB  "•l**"  ••••-/✓/-  &Ii<»LALE*»CNT 


Finite  element  arrangement  Fig-  21  Yielding  sequence  of  the  edges  in  th 

Prandtl  punch;  smooth  punch  material 

Slip  mechanism  small  node  contact 
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